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Abstract 

We recalculate the proton Dirac form factor based on the perturbative 
QCD factorization theorem which includes Sudakov suppression. The evolu- 
tion scale of the proton wave functions and the infrared cutoffs for the Su- 
dakov resummation are carefully chosen, such that the soft divergences from 
large coupling constants are diminished and perturbative QCD predictions 
are stablized. We find that the King-Sachrajda model for the proton wave 
function leads to results which are in better agreement with experimental 
data compared to the Chernyak-Zhitnitsky wave function. 
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1. Introduction 



Since the proposal of the improved perturbative QCD (PQCD) factoriza- 
tion formulas for exclusive processes, with the Sudakov resummation taken 
into account there have been many applications in the literature, such 
as the pion form factor 0, photon annihilation into pions the proton 
form factors |^, pion Compton scattering 0, proton-anti-proton annihi- 
lation 1^, and proton-proton Landshoff scattering These studies show 
that in the pion case the nonperturbative contributions from the end points 
of parton momentum fractions are moderated by Sudakov suppression, and 
perturbative predictions become relatively reliable. However, in the processes 
involving protons, because more partons share the external momentum, the 
infrared divergences associated with soft partons, which appear in hard scat- 
tering subamphtudes, are severer. It is then a concern whether Sudakov sup- 
pression of the end-point nonperturbative enhancements is strong enough to 
maintain the applicability of PQCD to the proton form factor at currently 
accessible energy scales. 

The improved factorization formalism has been applied to the proton form 
factor 1^ . However, the choice of the infrared cutoffs for the resummation was 
criticized [Q: The end-point enhancements are in fact not diminished com- 
pletely by Sudakov suppression under the above choice of cutoffs, implying 
that PQCD predictions remain unreliable. A modified choice of the cutoffs 
has been proposed 0, and the soft enhancements were found to be sup- 
pressed. Unfortunately, it turned out that the PQCD contributions amount 
only to half of the data, and hence it was concluded that higher-order or 
higher- twist corrections may be important 0. 

In this letter we shall recalculate the proton Dirac form factor based on 
the work of ^ by slightly modifying the infrared cutoffs for the resummation, 
and employing the more complete two-loop expression of the Sudakov factor. 
It will be shown that the end-point sensitivity is removed, and the PQCD 
predictions from one of the currently available models of the proton wave 
function match the experimental data well. We then confirm the applicability 
of the improved PQCD formalism for momentum transfer around few GeV. 
However, we emphasize that the uncertainty involved in our analysis is not 
negligibly small, and that the method in based on the overlap integral of 
the proton wave functions may be regarded as a complementary approach to 
ours. 
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2. Factorization 



According to the PQCD theory for exclusive processes |T^, the proton 



Dirac form factor, can be factorized into two types of subprocesses: wave 
functions which contain the nonperturbative information of the initial- and 
final-state protons, and a hard subamplitude which describes the scattering 
of a valence quark of the proton off the energetic photon. The former can 
not be calculated perturbatively, and needs to be parametrized by a model 
or to be derived by nonperturbative methods such as QCD sum rules. The 
latter, characterized by a large momentum flow, is calculable in perturbation 
theory. We quote directly the factorization formula for the proton form factor 
derived in 0: 

•J 

xHa'pi^'ap^{ki,k[,Q,^)Yo,p^{ki,P,n) , (1) 

with 

3 

{dx) = dxidx2dx2,5(^^Xi — 1) , 

1=1 

3 

{dkr) = dkiTdk2Tdk3TSC^kiT) . (2) 

1=1 

P = (P"'",0,0) is the initial-state proton momentum, and Xi = kf /P^ and 
kiT are the longitudinal momentum fraction and transverse momenta of the 
parton i, respectively. The primed variables P' = (0,P~,0), x[ = k~/P~ 
and k'-rp are associated with the final-state proton. = 2P ■ P' is the 
momumtum transfer. In the Breit frame we have P^ = P'^ = Q/\^. The 
scale fi is the renormalization and factorization scale. 

The initial distribution amplitude Yajs-y, defined by the matrix element of 
three local operators in axial gauge [0, |T^ , is given by 



= ^^[{VCU{l5N),V{h,P,f,) + {f^,CUN,A{h,P,f,) 

- {crlp''CU{ri^N),T{h,P,ii)] , (3) 
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where A''^ = 3 is the number of colors, \P) the initial proton state, u and d the 
quark fields, a, b and c the color indices, and a, /? and 7 the spinor indices. 
In our notation, 1,2 label the two w-quarks and 3 labels the ci-quark. The 



second form shows the explicit Dirac matrix structure [^, where /at is the 
normalization constant [|l^], the proton spinor, C the charge conjugation 
matrix and a^i, = [7^,7j,]/2. The amplitude Ya'/s'Y^k'i, P' , fi) for the final- 



state proton is defined similarly. By using the permutation symmetry |11 



and the constraint that the total isospin of the three quarks is equal to 1/2, 
it can be shown that the three functions V, A, and T are not independent, 
and related to a single function ip by [|ri[ 

1 

A{ki, ^2, /C3, P,fi) = - [ijj{k2, ki, fcs, P, n) - ij{ki, k2, fcs, P, /^)] , 

T{ki, k2, ks, P,f^) = \ [^{ki, h, k2, P, fi) + Hh, ks, fci, P, fi)] . (4) 

The hard subamplitude Ha'p'^'ap-y is obtained from the photon-quark scat- 
tering diagrams, and the expressions for the integrands Yaip'-yiHa'pi^iap-yYap-y 
are referred to Table I in Employing a series of permutations of the 
parton kinematic variables, Eq. (|l]) in Fourier transform space reduces to 

FliQ") = j:^ j\dx){dx'){dh)[fMr 

xHj{xi, x[, hi, Q, hi, P, P', /i) , (5) 

with bj the conjugate variable to kj^- and (dh) = dhidh2/ (27r)^. The explicit 
expressions of Hj and of \E'j in terms of ip will be below. 



3. Sudakov Suppression 

The Sudakov resummation of the large logarithms in ip leads to 



ilj{xi,hi,P,fi) = exp 



s{xi, w,Q)-3 / — 7g 



X(l){Xi,w) 



(6) 
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where the quark anomalous dimension 75(0:5) = —a. si tt in axial gauge governs 
the renormalization-group (RG) evolution of t/'. The function 0, obtained by 
factoring the Q dependence out of corresponds to the standard parton 
model. The exponent s is written as lir 



s(x, w, Q) 



■^Ql^ dp 
P 

where the anomalous dimensions A to two loops and B to one loop are 



In (^^j A{aM) + B{as{p)) 



(7) 



A 
B 




/3oln 



(8) 



Hf = 3 being the number of flavors, and 'Je the Euler constant. The two-loop 
running coupling constant, 



1 



Pi lnln(^VA2^ 



TT 



with the coefficients 



/?oln(^VA2) /3i ln2(/iVA2) ' 
33 -2nf 153 - 19nf 
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Pi 



24 



(9) 



(10) 



and the QCD scale A = Aqcd, will be substituted into Eq. (|^). 

The infrared cutoff w is chosen to be the inverse of a typical transverse 
distance among the three valence quarks. We try different definitions of this 
cutoff to determine its infiuence on the final result. One possible choice is 
w = l/bjnax, bmax = max(6/), / = 1,2,3, adopted in [|], with 63 = |bi - b2|. 
As long as all of these mass scales are much larger than A, the Sudakov 
form factor should not give any suppression. As one of these scales gets 
close to A, the Sudakov form factor tends to zero and suppresses this region. 
We find that choosing the infrared cutoff in this fashion suppresses all the 
infrared divergences and leads to a self-consistent calculation of the form fac- 
tor. However, this choice does not always correspond to a typical size of the 
three quark system. A more appropriate definition is obtained by considering 
it as a quark-diquark like configuration. The diquark constituents are taken 
to be those two quarks that are closest to each other in the transverse plane. 
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Figure 1: The typical transverse distance, dtyp.' The transverse distance be- 
tween the quarks A and B is the smallest among the three quarks. The diquark 
constituents are therfore considered to be the quarks A and B. The center of 
mass of the diquark, COM^q, is taken to be the central point of the line 
that connects these two quarks, dtyp is then defined as the distance between 
COMfiq and the third quark C. 



We now define the typical transverse distance, dtyp, as the distance between 
the center of mass of the diquark and the remaining third quark (Fig. p. 
This is clearly a more reasonable measure of the distance in the three quark 
system that can be resolved by a gluon. We shall therefore take the infrared 
cutoff as cw, where c is a parameter which is allowed to deviate slightly from 
unity. When we put c = 1, we recover the original choice of the cutoff. The 
introduction of this parameter c is natural from the viewpoint of the resum- 



mation, since the scale cw, with c of order unity, is as appropriate as w 113 



We choose c such that for a large number of randomly chosen triangles, of 
the type shown in Fig. |l|, we get for the average (dtyp/bmax) = l/c Defining 
c in such a way, gives c ~ 1.14. 

We find that both of these choices of the cutoff, lead to self-consistent 
calculations of the form factor in the sense that the form factor saturates at 
the large distance cutoff be- Remarkably, we find that with the small mod- 
ification of w into cw, which differs from what was used in 0, the results 
are in good agreement with experimental data. The dependence of the final 
results on the precise choice of scale cw shows that large distance contribu- 
tions cannot be completely dismissed, and give about 25-50% contribution 
at laboratory energies. Nevertheless, we find it encouraging that a physically 
motivated cutoff gives good agreement with experiments. 

The choice of scales for the Sudakov resummation in Eq. (0) is compared 
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to that adopted in |Q, where the different cutoffs bi are assigned to each 
exponent s and to each integral involving 7^: 



ilj{xi,hi,P,fi) = exp 



s{xi, hi, Q) + / — 7g 
1=1 \ -'^i 1^ 



X(l){Xi,w) . 



The Sudakov factor in Eq. (0) does not suppress the soft divergences from 
bi —>■ 1/A completely. For example, the divergences from bi — *• 1/A, which 
appear in (f){xi,w) at w A, survive as Xi 0, since s{xi,bi,Q) vanishes 
and s{x2, 62, Q) and s(x3, 63, Q) remain finite. On the other hand, w should 
play the role of the factorization scale, above which QCD corrections give the 
perturbative evolution of the wave function ip in Eq. (P), and below which 
QCD corrections are absorbed into the initial condition 0. It is then not 
reasonable to choose the cutoffs bi for the Sudakov resummation different 
from w. 



4. RG Evolution 

The RG evolution of the hard scattering subamplitudes is written as 

Hj{xi,x[,hi,Q,n) = exp -3^ / — 7g(a,(/i)) 

1=1 "'A' 

xHj{xi,Xi,hi,Q,tji,tj2) , (12) 
where the explicit expressions of t are 



;i-xi)(i-x;)Q,i/6i 



til = max 

^21 = max ^xix'iQ, 1/bi 
ti2 = t22 = max 



X2X'2Q, I/&2 



(13) 



The first scales in the brackets are associated with the longitudinal momenta 
of the exchanged gluons and the second scales with the transverse momenta. 
The arguments tji and tj2 of Hj denote that each is evaluated at the 
largest mass scale of the corresponding gluon. 
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Inserting Eqs. and (0) into Eq. (|^), we obtain 

~[ 27 JO JO JO 

xHj{xi, x'i, bi, Q, tji, tj2) ^j{xi, cw) 

X exp [-5'(xi,x'i,ciu,(5,tji,tj2)] , (14) 
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i^i = 3«.(tii)«s(ti2)^o ( V(l - 3^1) (1 - ^'i)Q^i ) ^0 ( \/x2X'2Qh ) 



with 



= 2as(i2i)tts(i22)^o (^Va:ix'iQ6ij i^o (^Va;2a;2Q^2j • (15) 

The variable 6 is the angle between bi and b2. Kq is the modified Bessel 
function of order zero. The expressions for are 



^1 



2(00')l23 + 8(rT')l23 + 2(00')l32 + 8(rT')l32 " (00')321 " (00')231 



'l-xi){l-x[) 



_ 2(00Oi32-2(TTOi23 (00Ol23-8(TTOi32- 2(000321 . . 

(l-X2)(l-x'0 (i_a;3)(l-xi) ' ^ ' 

which group together the products of the initial and final wave functions in 
the notation 

(00')i23 = 0(a;i, X2, X3, cw)(f){x[, X2, Xg, cw) . (17) 

[TT') is defined similarly based on Eq. (^ but with ip replaced by 0. The 
Sudakov exponent S is given by 

S{xi,x'i,cw,Q,tji,tj2) = '^s{xi,cw,Q) + 3 —'jqiasifi)) 

Jew fi 



j2s{x'i,cw,Q)+3 r^"^7ja,(^)) .(18) 

l—-^ Jew H 



For the wave function 0, we will consider both the Chernyak-Zhitnitsky 
(CZ) model and King-Sachrajda (KS) model [|I4|. They are decomposed 
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in terms of the first six Appel polynomials Aj{xi), which are eigensolutions 



of the evolution equation for the nucleon wave function |jTO|, [15 



{Xi,w) = (l)as{Xi)J2^j 
j=0 



"^(/Lto) 



(19) 



with /io ~ 1 GeV. The constants Nj, aj and bj are given in Table I. (pasi^i) = 
120x1X2X3 is the asymptotic form of 0. The evolution of the dimensional 
constant /at is given by 



fNiw) = /Ar(/io) 



l/(6/3o) 



(20) 



with f^ifio) = (5.2 ± 0.3) X 10-3 GeV^ |in|. 



5. Numerical Results 



In order to be able to calculate the seven-dimensional integral, Eq. (14), 
we use the VEGAS Monte Carlo routine fl^. We set the Sudakov factor 
exp(— 5*) to unity in the small b region where it displays a small enhance- 
ment, since in this region higher-order corrections should be absorbed into the 
hard scattering instead of into the wave function, giving its evolution. 
Therefore we have also set the factor exp[—s{C,,cw,Q)] to unity whenever 
C,Q/\/2 < cw. As cw approaches A, the Sudakov factor vanishes, implying 
that the whole integrand of Eq. (|14D also vanishes. 

First we choose the parameter value c = 1. The results of Q^Ff for 
A = 0.2 GeV from the use of the KS wave function, along with the experi- 
mental data 0, |l8l, are shown in Fig. The PQCD predictions amount 
only to about 60% of the data. It is then possible that higher-order or higher- 
Fock-state contributions are necessary for the explanation of the data, which 
are certainly worth of further studies. However, before jumping to that con- 
clusion, we investigate the effect from the freedom of varying the parameter 
c. The results with c = 1.14 are also displayed in Fig. |^. It is found that 
our predictions match the data well. Note that varying c makes a difference 
in the resummation at the level of next-to-leading logarithms, which can be 
regarded as an uncertainty of our formalism. Therefore, we argue that the 
current data can be explained within the uncertainty of our approach. 



9 



2 




Q I \ \ \ \ \ \ \ I 

5 10 15 20 25 30 35 40 
q2 (GeV2) 

Figure 2: Dependence of Q^Ff on for the use of the KS wave function 
(c = 1.14, solid line; c = 1, dense-dot line) and for the CZ wave function 
(c = 1.14, dashed line; c = 1, dotted line). The experimental data with 
errorbars are also shown. 

Following 1^, we should analyze how the contributions to Q'^Ff are dis- 
tributed in the &1-62 plane. The integration is done with both variables bi 
and 62 cut off at a common value be- If the perturbative region dominates, 
most of the contributions will be quickly accumulated below a small be. The 
numerical outcomes (with c = 1.14) are shown in Fig. ^. All the curves, 
showing the dependence of Q^Fi on b^ increase from the origin and reach 
their full height at b^. = 0.9/A. The curves exhibit small humps at the high 
end of be, which imply that the evolution of the wave function gives a small 
negative contribution in the large b region. A standard of self-consistency is 
that 50% of the whole amount of Q'^Ff is accumulated from the region with 
a^/TT smaller than 0.5. Based on this standard, the results with > 10 
GeV^ are reliable. Therefore, the applicability of PQCD to the proton form 
factor at currently accessible energy scale ~ 35 GeV^ is justified. 
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bcA 

Figure 3: Dependence of Q^Ff on the cutoff be with the KS wave function 
for Q^ = 12 GeV^ (dotted line), = iq QeV^ (dashed Une), = 2h GeV^ 
(dense-dot line), and = 36 GeV^ (solid line). 

The CZ wave function is also employed, and the corresponding results are 
shown in Fig. ^. It is observed that the values are only about 2/3 and 3/4 
of those derived from the KS model with c=1.14 and c=l respectively, and 
are far below the data. Hence, we claim that the KS proton wave function 
is more phenomenologically appropriate. 

6. Conclusion 

In this work we have modified the choice of the infrared cutoffs for the 
resummation, and employed the more complete two-loop expression of the 
Sudakov factor compared to the previous analyses. With these modifications, 
we have been able to explain self-consistently the experimental data of the 
proton Dirac form factor for > 10 GeV^, within the uncertainty of our 
formalism. We should emphasize that though the nonperturbative region 
denoted by 6 — 1/A does become less important in our analysis, the coupling 
constant is not so small that we could consider the perturbative results 
as exact. Therefore, nonperturbative contributions may be comparable to 
the perturbative ones at the currently accessible energies. A complementary 
study based on nonperturbative approaches such as QCD sum rules and 
the determination of the transition of the proton form factor to PQCD, as 
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performed in [0, are then essential. 

The analysis presented here is not conclusive even in the PQCD frame- 
work. The uncertainty at the level of next-to-leading logarithms indicates 
that higher-order corrections to the evolutions of the wave function and of 
the hard subamplitudes need to be evaluated. The contributions from higher 
Fock states should be investigated, which may be important in the interme- 
diate energy range. 

It is found that the KS wave function is more phenomenologically appro- 
priate, which will be adopted in the future studies of QCD processes involving 
protons. 
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Table I. Appel polynomial coefficients in Eq. 
function (f){xi,w) of the CZ and KS models [ pi] , 
GeV 0. 



191) for the nucleon wave 
^ with the scale /xq ~ 1 



J a,(CZ) a,(KS) N, b, A,(x.) 

1.00 1.00 1 1 

1 0.410 0.310 21/2 20/9 xi - X3 

2 -0.550 -0.370 7/2 24/9 2 - 3(xi + a;3) 

3 0.357 0.630 63/10 32/9 2 - 7(xi + X3) + 8(x? + x^) + 4x1X3 

4 -0.0122 0.00333 567/2 40/9 xi - X3 - (4/3)(x? - xi) 

5 0.00106 0.0632 81/5 42/9 2 - 7(xi + X3) + 14xiX3 
+(14/3)(x^ + xi) 
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